Ferromagnetic Relaxation by Magnon Induced Currents 
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A theory for calculating spin wave relaxation times based on the magnon-electron interaction is 
developed. The theory incorporates a thin film geometry and is valid for a large range of magnon 
frequencies and wave vectors. For high conductivity metals such as permalloy, the wave vector 
dependent damping constant approaches values as high as 0.2, showing the large magnitude of the 
effect, and can dominate experimentally observed relaxation. 



One of the fundamental problems of magnetism is to 
determine a mechanism for dissipation of energy in a 
system subject to a change in the direction of the ex- 
ternal magnetic field. Historically, these ferromagnetic 
relaxation processes have been explored by ferromag- 
netic resonance (FMR) in which the absorption of a small 
RF frequency field applied perpendicular to a large DC 
field is measured. More recently, direct measurement of 
large angle switching has been made by multiple groups 
jH EE El • Time resolved Kerr microscopy has made it pos- 
sible to track isolated magnetic relaxation processes with 
picosecond temporal resolution Q] and also to image the 
individual components of the precessing magnetization 
3J. It has been previously shown that for many exam- 
ples of large angle switching, particularly those involv- 
ing materials with large magnetizations, the dominant 
relaxation process is very different than that applicable 
to FMR. In particular, it was found that the coherent 
mode scatters with thermal magnons to form two k ^ 
magnons p|. This 4- magnon process rapidly escalates 
as the k ^ magnon levels are populated, thus pro- 
moting additional scattering. This previous work thus 
accounted for the rapid movement of the magnetization 
into the new direction, but the dissipation of the mag- 
netic energy stored in the k ^ modes still must be 
addressed. 

In a conducting ferromagnet the interaction between 
the conduction electrons and the magnons become very 
important. The magnetic field generated by the spin 
wave is time dependent and therefore, by Faraday's law, 
it creates an electric field in the system. These electric 
fields, unlike conventional eddy currents, are wave like in 
nature. In a metallic system, the fields drive the conduc- 
tion electrons. These magnon induced currents help dis- 
sipate the energy of the system by Joule heating. Abra- 
hams Q addressed this question half a century ago by 
taking into account the interaction between spin waves 
and conduction electrons. However, his bulk estimates 
predicted relaxation times one order of magnitude less 
than required to explain FMR line widths. Subsequently 
several attempts have been made to give a consistent the- 
ory of ferromagnetic relaxation from the point of view of 
an FMR experiment Q. Later, Almeida and Mills @ 
explored the same interaction and derived the Green's 



function for the limited case of small angle precession in 
the absence of quantum mechanical exchange, i.e. in the 
long wavelength limit. 

In the present work, we avoid much of the limitations 
found in previous work. Solution of the general prob- 
lem is difficult because the magnon-induced currents gen- 
erate new fields which further affect magnons and cre- 
ate new currents. Here we show that expansion in the 
small parameter 4iruja/c 2 k 2 allows an explicit solution to 
the general problem. It allows prediction of decay rates 
for magnons of arbitrary frequency and large amplitudes 
that are limited only by the Holstein Primakoff transfor- 
mation In particular, it can be applied to the problem 
of magnons generated by four-magnon scattering after a 
large angle rotation such is common in modern switch- 
ing experiments and technological applications such as 
magnetic recording. We also discuss our results in con- 
text of the spin wave resonance experiments capable of 
measuring the linewidth of the higher order k modes. 
Historically, these experiments were used to measure the 
exchange constant of the material. 

We consider an infinite film of thickness d made of 
ferromagnetic metal. The top and the bottom surfaces 
of the film are at z = d and z — respectively (see Fig. 
^) . A spin wave of wave vector k = kx and frequency u 
is excited in the system. 

We write the electric and magnetic fields in the system 
as a series expansion 
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where a is the conductivity of the medium and c is the 
velocity of light. For typical frequency and wave vector, 
this expansion parameter is quite small, e.g. 10~ 4 for 
Fe (<7 = 9 x 10 16 s _1 ). Therefore, the series converges 
rapidly and only the leading term has practical interest. 
The n— th order terms in the expansion of Eq. Q obey 
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Figure 1: Schematic diagram showing the geometry used in 
the paper. An infinite ferromagnetic metallic slab with thick- 
ness d along the z-direction is magnetized in-plane by apply- 
ing an external field H. We consider spin waves with wave 
vector k propagating parallel to the surfaces. 



the Maxwell equations 

v-eW = 

V • (H (n) + 47rM<5 n0 ) = 



V x E (n) = - 



1 dB^ 

c dt 



VxH (n+l) = l^ E (n) 



(2) 



We neglect the displacement current in the last expres- 
sion owing to oj <C a. From Eq. J5J V x H' ) = 0, so we 
can write the magnetic field as the gradient of a magnetic 
scalar potential: H' ) = — V</>m- The scalar potential has 
a volume and a surface term and the zeroth order mag- 
netic field produced by the spin waves can be written as 



H (°) = - d 3 r 



,V-M(r')(r-r') 



iV n-M(r')(r-rQ 



(3) 



where M is the magnetization of the sample and n is the 
outward normal to the surface carrying magnetic charge. 
We shall consider films to be thin enough that the spin 
waves are confined to the x — y plane only. We shall 
consider two specific cases: k_LM and k||M. 

We consider the spin wave propagating along the x 
direction in the thin film. In configuration I, the magne- 
tization is precessing in the x — z plane 



M 



Moy + e [cos (kx — tot) x + sin (kx — tot) z] 



where Mq is the component of magnetization perpendic- 
ular to the plane of precession and e is the amplitude of 



precession. In our calculation, we need not restrict e to 
be small compared to Mq. The zeroth order magnetic 
field from Eq. © for this configuration is 



Hl {0) = 4e7T cos(kx - ut) 
= 4e7r cos(kx — ujt) 
Hf> = 

#1(0) = „ 4e?r S [ n (kx - LUt)t 
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According to Eq. J5J H 1 ^ ) generates an electric field 
E 1 ^ ) which has only one nonzero component 
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where £ = 2neui/c. Note the asymmetry of the solution 
for positive and negative values of k. The profile of the 
electric field for k > and k < are mirror symmetric 
with respect to z = d/2. 

In configuration II, the magnetization is precessing in 
the y — z plane 

M 11 = Mqx + e [cos (kx — ut) y + sin (kx — tot) z] 

Since V • M n = only the surface term of Eq. © con- 
tributes to the magnetic field given by (for k > 0) 



-2e7rcos(fcx — ujt)[e 
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H?<® = -2ensm(kx-u;t)[e kiz - d) + e~ kz ] 
which generates an electric field 
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where A = (£/fc)[(l - e~ 
(£,/k)\(e kd — 1)/ sinh(fcrf)]. Note that unlike in the previ- 
ous configuration the components of the electric field are 
symmetric with respect to the two surfaces. 

The energy stored in the form of spin waves is dissi- 
pated from the system by the current generated by the 
electric field induced by the precessing spins. This in- 
duced electric field drives the free electrons in the metal 
to produce the magnon induced current. The ohmic 
power loss per unit volume due to this magnon induced 
current can be written as 



P = lim - ! 
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Figure 2: Dispersion relation for configurations I and II ac- 
cording to Eqs. JHJ and ® for an infinite iron slab using 
M s = 1700 emu/cc, A = 2 x 10" 6 erg/cm, 7 = 1.76 x 10 7 
(Oe s)" 1 and H eKt = 1000 Oe. 



Integrating over the square of the electric field described 
in Eqs. (Jlj and J^}, we obtain the power dissipation per 
unit volume 



P\k,d) 
P n (k,d) 



= 2/3 [-3 + 2dk + 4e- dk - e- 2dk ] 

a„-2kd 

= pe _ [(— 15 + 8dk)e 3dk 

+ (9 + \Qdk)e 2dk + (7 + 2dk)e dk - 1] 



where 



2ov 2 e 2 u 2 
dk 3 c 2 



(6) 



(7) 



The power dissipation clearly depends on u> which is 
a function of k. Owing to the strong influence of the 
magnetostatic energy within the magnon Hamiltonian, 
the derivation of this relationship is nontrivial, but, for- 
tunately has been described by previous workers. Es- 
sentially, the Hamiltonian of the system has contribution 
from exchange, magnetostatic and Zeeman energy. We 
shall restrict ourselves to isotropic systems and there- 
fore crystallographic anisotropy will have negligible ef- 
fect to our result. We will only consider magnons with 
wavelength much greater than the lattice constant, e.g. 
k < 10 7 cm -1 which applies to most magnons of interest. 
The dispersion relation for a thin film is given by [lH 



uj = 7 
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where Hi = H cxt — NM$ is the internal field, A is the 
exchange constant, N is the demagnetizing factor and 
-ffext is the external magnetic field. The magnetostatic 
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Figure 3: Relaxation time for an infinite iron thin film with 
thickness di = 1 X 10~ 6 , d 2 = 5 x 10" 6 , d 3 = 1 x 10" 5 and 
c?4 = 5 x 10~ 5 cm for (a) configuration I and (b) configuration 
II. 



contribution is given by 
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The dispersion relation for an iron thin film of various 
thicknesses is shown in Fig. [21 Note that the curves 
converge for k Si 2 x 10 6 cm -1 where exchange interaction 
starts dominating the thickness dependent magnetostatic 
interaction. 

The energy density of the system consists of the ex- 
change, the magnetostatic and the Zeeman terms 



£(k,d) = A 



fVM 



£ d - M • H e 



(10) 



where the magnetostatic contribution to the energy is ob- 
tained by taking the negative scalar product of magneti- 
zation and the magnetic field produced by the magnons 
and is given by 



2we 2 
2tt£ 2 
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Eqs. ©, Q) ©, ©, JTUJl and (JTTI show that the energy 
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Figure 4: The quantity 1 /tuj as a function of the wave vector 
for two different thicknesses. 



decays exponentially with relaxation time: 



r(k,d) 



S(k,d) 
P(k,d) 



This is the leading order term contributing to the relax- 
ation time as the energy of the system gets renormalized 
by the magnetic field generated by the magnon-induccd 
currents. Fig. |21 shows the relaxation time of spin waves 
in an iron thin film where the jSpin wave is confined in 
the x — y plane. According to [4j the four magnon scat- 
tering produces magnons with wave vectors in the range 
1 — 5 x 10 6 cm -1 , ft is important to note that magneto- 
static is the dominant interaction in the long wavelength 
limit whereas exchange takes over in the short wavelength 
regime. The crossover, which happens around k ~ 10 6 
cm -1 in the range of thickness of the sample we are in- 
terested in, is of the same order of magnitude for bulk 
magnetic crystals 0. This feature is nicely reproduced 
in our result where changing the value of the exchange 
constant only affects the curves for k > 10 6 cm -1 . 

Both for configurations I and II, t increases with k for 
thicker films in the magnetostatic regime and saturates in 
the exchange regime. For such films the electric field and 
hence the power shows an inverse dependence with the 
wave vector in the long wavelength limit whereas the en- 
ergy of the system is independent of k in configuration I 
and changes slowly with k in configuration II. This makes 
r an increasing function of k except for small thicknesses 
of the film in which case the power dissipation becomes 
nearly independent of k. However, both £ and P in- 
creases as k 2 in the exchange regime thereby making r 
independent of k. 

Fig. 01 shows the behavior of the quantity 1/rw as a 
function of the wave vector. This quantity, which can 
be interpreted as a wave vector dependent damping con- 
stant, is a measure of the amount of energy taken away 
from the system per precessional cycle. This damping 
constant should not be confused with the typical Gilbert 



damping constant used for uniform rotation of the el- 
ement magnetization. However, the presence of values 
approaching 0.2 illustrates the large magnitude of the ef- 
fects described here. 

Observation of spin wave resonance 13] in a thin film 
allows calculation of relaxation times from measured line 
widths. For example, analysis of Okochi's 0] data for his 
first spin wave resonance mode yields a damping constant 
a = 0.0063. Calculations for this configuration yield P = 
32-rT 2 e 2 Lu 2 a/(c 2 k 2 ) and£ = (2Ak 2 /Mj-4TT+H cxt /M s )e 2 . 
The resulting value of the damping constant is 0.0095 
(using A = 10~ 6 erg/cm and a = 2.9 x 10 16 s" 1 [3 for 
FeNi) . The discrepancy is presumably due to conductiv- 
ity differences between the two thin film samples. It is 
worth noting these experiments rely on exciting stand- 
ing waves along the perpendicular direction of the thin 
film. This geometry typically minimizes the effect (rela- 
tive to the configurations discussed elsewhere in this text) 
of magnon induced currents in the relaxation time be- 
cause the magnon generated magnetic field is zero which 
makes the induced electric fields weaker and the rate of 
energy dissipation slower. 

The effect of conductivity on the magnon-electron dis- 
sipation mechanism can be studied in magnetic semicon- 
ductors such as CdCr2Se4 and HgCr2Se4 whose conduc- 
tivity can be tuned by the amount of Ag doping |l6t Il7| . 
The typical conductivity of such materials is several or- 
ders of magnitude lower than that of a ferromagnetic 
metal. For example, a 0.75 mole % Ag doped CdCr2Se4 
has a = 4.5 x 10 11 s" 1 at T = 120K. Therefore, the co- 
efficient of expansion in Eq. Q becomes much smaller 
compared to that for a ferromagnetic metal thereby mak- 
ing our formalism highly applicable for such materials. 
Assuming a resonance field H Tes = 3500 Oe 0] and 
lu = 7-f/j-cs we obtain a linewidth AH =1/ (tj) = 78 
Oe for configuration I and AH — 144 Oe for configura- 
tion II, using d = 0.2 mm and k = l/d. These values 
are in very good agreement with the FMR linewidths ob- 
served by Ferreira and Coutinho-Filho Under the 
experimental condition the exchange contribution to the 
energy is negligible and the magnetostatic and the Zee- 
man energies are of the same order of magnitude. In 
this limit, therefore, the magnon-electron contribution to 
the energy dissipation (calculated here) is comparable to 
that of conventional Eddy current loss obtained from the 
FMR linewidth by subtracting the effect of two magnon 
scattering [lij ]. 

We conclude by proposing the following picture of fer- 
romagnetic relaxation in switching experiments. We ex- 
pect that the initial rapid approach of magnetization 
direction to equilibrium is enabled by magnon-magnon 
scattering that converts the energy into the higher spin 
wave modes. These modes then decay at a slower pace via 
the magnon-electron interaction described here or by the 
traditionally invoked mechanisms in less pure, lower con- 
ductivity films. This delay will lead to a small reduction 



in magnetization which appears to have been observed 
by Silva et al [Tg| . 
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